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SUMMARY 


A study has been made of the factors which place limits on the theoretical 
leading-edge thrust, and an empirical method for estimating attainable thrust 
has been developed. The method is based on the use of simple sweep theory to 
permit a two-dimensional analysis, the use of theoretical airfoil programs to 
define thrust dependence on local geometric characteristics, and the examination 
of experimental two-dimensional airfoil data to define limitations imposed by 
local Mach numbers and Reynolds numbers. The applicability of the method was 
demonstrated by comparisons of theoretical and experimental aerodynamic charac- 
teristics for a series of wing-body configurations. Generally, good predictions 
of the attainable thrust and its influence on lift and drag characteristics were 
obtained over a range of Mach numbers from 0.24 to 1.3. 

The method is compatible with the Polhamus leading-edge suction analogy for 
fully detected vortex flow, when the analogy is considered to represent the 
limiting condition of a gradual rotation of the suction vector as leading-edge 
thrust is lost. An additional advantage of the method is the possibility it 
presents for designing wings in an iterative manner to maximize thrust and the 
attendant performance benefits. 


INTRODUCTION 

Leading-edge thrust is an important but little understood aerodynamic phe- 
nomena that can have a large influence on wing aerodynamic performance. This 
force results from the high local velocities and the accompanying low pressures 
which occur as air flows from a stagnation point on the undersurface of the wing 
around the leading edge to the upper surface. At low subsonic speeds, the high 
aerodynamic efficiency of uncambered wings with high aspect ratios depends 
directly on the presence of leading-edge thrust to counteract the drag arising 
from pressure forces acting on the remainder of the airfoil. Leading-edge 
thrust tends to diminish with increasing speeds, but may be present to some 
degree even in the supersonic speed regime, provided the leading edge is swept 
behind the Mach angle. 

Leading-edge thrust for subsonic speeds may be predicted by a variety of 
methods including a vortex-lattice computer program (ref. 1) capable of handling 
wings of complex planform. At supersonic speeds, leading-edge thrust for wings 
with straight leading edges may be determined by purely analytic means (e.g., 
ref. 2). A recently developed computer method (ref. 3) has extended this capa- 
bility to wings of arbitrary planform with twist and camber. These methods, 
however, provide estimates of only the theoretical thrust which may or may not 
be attainable in the real flow. 

At present, methods for estimating the fraction of the thrusting force 
actually attainable (e.g., ref. 4) are based on average conditions on the wing 
as a whole and are applicable at or near cruise conditions. Thus, the existing 



methods provide no information on locally attainable thrust or on its spanwise 
distribution. Such information would be useful in understanding the phenomena 
and could lead to the development of design methods for maximization of attain- 
able thrust and the attendant performance benefits. 

This report presents a study of the factors which place limits on the 
thrust and describes an empirical method for estimation of attainable thrust. 
The method is suitable for programming as a subroutine in existing methods for 
estimation of the theoretical leading-edge thrust. 
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SYMBOLS 


wing span 

local wing chord 

mean aerodynamic chord 

average wing chord, S/b 

theoretical section thrust coefficient, dt/dy qc 
attainable section thrust coefficient, dt*/dy qc 
axial or chord force coefficient 

increment of axial force coefficient due to leading-edge thrust 

drag coefficient 

drag coefficient at zero lift 

lift coefficient 

normal force coefficient 

increment of normal force coefficient due to rotation of leading-edge 
suction vector 

pressure coefficient 

limiting pressure coefficient used in definition of attainable thrust 
vacuum pressure coefficient, -2/yM^ 


theoretical wing thrust coefficient. 


attainable wing thrust coefficient. 
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k-L,k 2 ,k 3 ,k 4 constants used in airfoil section definition 

Krp fraction of theoretical thrust actually attainable , c */ c t; 

l wing overall length 

L/D lift-drag ratio 

M free- stream Mach number 

M e equivalent Mach number used in definition of %« 

q dynamic pressure 

r leading-edge radius 

R free-stream Reynolds number based on c 

s distance along wing leading edge 

5 wing area 

t theoretical section leading-edge thrust 

t* attainable section leading-edge thrust 

x,y,z Cartesian coordinate system 

a angle of attack , radians unless otherwise specified 

y ratio of specific heats 

6 angle between tangent to local camber surface and wing-chord plane, 

positive for trailing edge up, deg 

V\ location of maximum wing section thickness as fraction of chord 

A le leading-edge sweep angle, deg 

A m maximum- thickness-line sweep angle, deg 

A te trailing-edge sweep angle, deg 

T maximum wing section thickness 

Subscripts : 

i - 1 , 2, 3, 4, . . . 

n quantities pertaining to wing section normal to leading edge with 

maximum thickness at midchord (see fig. 1) 
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METHOD DEVELOPMENT 


Development of this method for the prediction of attainable leading-edge 
thrust is based on the fundamental principle that such forces must result from 
pressures acting on a surface and limitations on the attainable pressures and 
the surface areas on which they act form the only restraints on the achievement 
of theoretical thrust. The general plan of development covers three phases. 
First, the relationships between streamwise airfoil sections and sections normal 
to the leading edge are established so that a two-dimensional analysis may be 
conducted. Then a program for subsonic airfoils is employed to define limita- 
tions on the theoretical thrust imposed by airfoil geometry restraints and by 
arbitrarily defined limiting pressures. Finally, experimental data for two- 
dimensional airfoils were used to evaluate the practical limiting pressures and 
their dependence on Mach number and Reynolds number. 


Normal Airfoil Derivation 

For purposes of this analysis, the familiar concepts of simple sweep theory 
were employed. As shown in figure 1, the free-stream flow is separated into two 
parts, with one component of velocity parallel to the leading edge and the other 
perpendicular or normal to the leading edge. It is assumed that, through the 
transformations to be discussed, the streamwise section pressure distribution 
and the leading-edge thrust characteristics can be related to the characteristics 
of a two-dimensional section which is normal to the leading edge and operates in 
a velocity field defined by the normal flow vector. To define the airfoil sec- 
tion normal to the leading edge at any given span station, a superimposed arrow 
wing is introduced. This phantom wing has the same sweep angles of the leading 
edge, the trailing edge, and the maximum-thickness line as does the actual wing 
at the same span station. Derivations of the following relationships between 
normal and streamwise quantities are given in the appendix. 

The normal-flow Mach number is defined as 


M n - M cos Aj_ e 


( 1 ) 


The chord of the normal section is defined so as to place the maximum thickness 
at midchord. This is done to reduce the number of variables which must be 
studied in a subsequent analysis of two-dimensional section pressure and thrust 
characteristics. The ratio of the normal section chord to the streamwise section 
chord may be expressed as 


c n _ 2r\ 

c sin A le Ql - n) tan + n tan A fc J + cos A le 


( 2 ) 


4 



For an unswept trailing edge, equation (2) simplifies to 
C n 2 n cos A le 

c i „ .2a 

1 - n sin A le 

And for unswept leading and trailing edges, or two-dimensional sections, the 
equation may be further reduced to 



The thickness-to-chord ratio of the normal section may be expressed as 

^n_ _ 1 

c n c 2rj cos A le 

and the ratio of leading-edge radius to chord for the normal section is 


^n = r 1 

c n c 2q cos 2 A le 


(4) 


The normal section thrust coefficient is related to the streamwise section 
thrust coefficient by 


c 


t,n 


c t 


_c 1 

c n cos 2 A^_ e 


The normal flow Reynolds number is 


n * 

R n = r -=- cos A 


le 


(5) 


( 6 ) 


In regions of the wing leading edge away from the apex and away from the 
wing-body juncture, the preceding expressions are believed to provide a reason- 
able basis for two-dimensional analysis of leading-edge thrust phenomena. In 
those regions where the analysis is most questionable, thrust values are gen- 
erally small and errors in the attainable levels should have little impact on 
the overall performance of the wing. 


Theoretical Airfoil Analysis 

For the series of symmetrical two-dimensional airfoil sections shown in 
figure 2, the subsonic airfoil program of reference 5 was employed to define 
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pressure distributions and integrated thrust coefficients. The airfoils were 
defined by the following equation: 

z = k^^fx + k 2 x + k^x 3//2 + k^x 2 


in which the coefficients were selected to produce the required chord, thick- 
ness, and leading-edge radius for a wing section with maximum thickness at mid- 
chord. Maximum airfoil thickness ranged from 6 percent to 18 percent of the 
chord, and leading-edge radius ranged from 2 percent to 16 percent of the maxi- 
mum thickness. For a given airfoil, pressure distributions and thrust coeffi- 
cients covering a range of normal Mach numbers were found by applying the 
Prandtl-Glauert rule 

Cp^l - M n ^ = Constant 


to pressure distributions obtained at 
a Mach number of 0.01. This simple 
means of handling Mach number effects 
was employed for the sake of consistency 
with methods using linearized theory for 
estimating theoretical thrust. A typi- 
cal pressure distribution for a large 
angle of attack is shown in sketch (a) . 
Program values of the integrated section 
thrust coefficient 

- £ f c p az 

were found to be relatively independent 
of the airfoil thickness and the leading- 
edge radius and to be in reasonably good 
agreement with the two-dimensional theo- 
retical value, 2ira^/\[l - M n 2 . The 
suction peak of the theoretical pressure 
distribution can be quite large, often 
exceeding the vacuum pressure limit for 
a given Mach number. Thus the theoreti- 
cal section thrust coefficient can be 
unrealistically high. 


z 




Sketch (a) 
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In order to determine the effect on 
the thrust coefficient of realistically 
attainable pressure distributions , the 
program integration was performed twice : 
once without limitation as shown in 
sketch (a) and once with limitation to 
values greater than or equal to arbi- 
trarily defined pressure coefficients 
C p,lim (limiting the suction peak) as 
shown in sketch (b) . This pressure 
limitation is intended to account, in 
an approximate way, for two of the 
factors which limit attainable thrust: 
the failure to attain theoretical peak 
suction pressures and the tendency of 
these peaks to occur at a more rearward 
position on the airfoil. Definition of 
values of the effective limiting pres- 
sure coefficient will be addressed in a 
later section of the paper. The limited 
value of c t is designated c£ and the 
attainable thrust ratio or the thrust 
factor is simply: K T = c*/c t . 

Shown in figure 3 is an example of the variation of the thrust factor with 
angle of attack (and with the theoretical thrust coefficient) for a given normal 
airfoil section at a given normal Mach number. Inset sketches show pressure 
distributions for 5°, 10°, 15°, and 20° angles of attack. In this example, the 
limiting pressure was set equal to the vacuum limit. As the figure shows, this 
limitation can be quite severe for large angles. 

Program data are shown in figure 4 for the complete range of airfoil param- 
eters and for normal Mach numbers of 0.3, 0.5, 0.7, and 0.9 with a limiting 
pressure equal to the vacuum limit. For a given normal Mach number, the thrust 
factor Krp was found to correlate well with the parameter 



A thrust factor which decreases with increasing theoretical thrust is clearly 
shown. The correlation plot (fig. 4) also shows the tendency of increased thick- 
ness or increased leading-edge radius to improve the thrust factor. The curve 
fit is described by a single equation covering the full range of normal airfoil 
parameters and normal Mach numbers : 



z 



z 


Sketch (b) 
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The limitation of K T to values no greater than 1.0 permits attainable thrust 
to equal but not exceed the theoretical thrust. Note that the curve fit allows 
the thrust to go to zero when either the thickness or the leading-edge radius 
goes to zero. If both thickness and leading-edge radius go to zero, zero thrust 
is to be expected. However, as long as there is some thickness, an airfoil with 
a leading-edge radius of zero could produce a small amount of thrust because of 
the upper-surface suction pressure peak acting on forward facing slopes. Prob- 
ably a more important consideration is that real leading-edge radii can never be 
zero and some attempt should be made to estimate an effective leading-edge 
radius for airfoils with theoretically sharp edges. 

Equation (7) was developed to 
cover data in which the limiting 
pressure coefficient was set equal 
to the vacuum pressure coefficient. 

The results, however, can be general- 
ized to cover a full range of limiting 
pressures between 0 and Cp, vac for a 
given normal Mach number M n by means 
of the following logic. As illustrated 
by the pressure distributions shown in 
sketch (c) for a given airfoil section 
at a given angle of attack, the pres- 
sure coefficient at any point on the 
airfoil will vary with Mach number 
according to the Prandtl-Glauert rule. P 
If the limiting pressure Cp 
changes in this same fashion, the 
thrust factor Kip will be the same at 
all Mach numbers. Thus for any Mach 
number M e different from the normal 
Mach number under consideration, Kip 
will be the same as for M n provided 
that 

Jl - M n 2 

Cp, lim = Cp,lim( M n) < — 

\l 1 " M e 2 

Then if M e is selected so that Cp is equal to the limiting vacuum 

pressure for that Mach number Cp vac (M e ) , the appropriate value of Kip for the 
normal Mach number under consideration can be found from equation (7) by substi- 


tuting M e for 

M n* 

The required M e 

can be found by setting c p /V ac^ M e^ 

equal to Cp flim 

(M e ) 

, the intersection 

point of the curves shown in sketch (c) , 

and solving for 

**e- 

Thus , 


c 

(M e ) 

” C p,lim^ M e^ 


p , vac 
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-2 


ll - M 


2 ^p,lim^n^ 


M 


and on solving for the equivalent normal Mach number. 


M e = 




^ C p, lim^ 


1 - M 


2 1 


^1 + ^YCp f ii m \|l M n ^ 1 


( 8 ) 


Thus, the thrust factor can be expressed as 


Kjp = 


2(l 



/r 

n f n 


c n \ c n 


0.4 


0.6 


M 


e 




but not greater than 1.0 


( 9 ) 


where M , as defined by equation (8), covers the full range of possible limit- 
ing pressures. These expressions (eqs. (8) and (9)) describe the variation of 
the thrust factor with the theoretical thrust, with airfoil geometric param- 
eters, with normal Mach number, and with limiting pressure coefficients, which 
will be defined in the next section. 


Experimental Airfoil Analysis 

In order to define practical values of the limiting pressure coefficient, 
the as yet incomplete prediction method was applied to experimental two- 
dimensional airfoil data (refs. 6 and 7) for symmetrical sections. Correlations 
of axial force coefficients predicted by this method with experimentally deter- 
mined axial force coefficients as shown in figures 5, 6, and 7 were used to 
determine, by iteration, values of c p f Hm appeared to match the experi- 

mental trends . 

Figure 5 is intended to show the leading-edge thrust behavior with changes 
in wing thickness for a low normal Mach number of 0.3. Although there are large 
changes in the fraction of theoretical leading-edge thrust actually achieved, 
the differences in the limiting pressure coefficient required for correlation 
are relatively minor, varying from 12.5 to 15.0 percent of the vacuum pressure 
coefficient C p /V ac for th i s Mach number. In order to show the sensitivity of 
the limiting pressure to uncertainties in the experimental data, curves of AC a 
as a function of a are shown in figure 5 for values of Cp f ]_±xa 20 P ercent 
greater than and 20 percent less than the chosen value. 
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Shown in figure 6 are the variations in leading-edge thrust with changes 
in the location of maximum thickness. These data have a somewhat higher level 
of indicated thrust than shown in figure 5, apparently caused by the use of a 
tunnel-wall bleed system having a large pressure equalizing duct rather than a 
small one. Note the differences in the data of figures 5 and 6 for the identi- 
cal NACA 64A009 airfoil section. Limiting pressure coefficients indicated by 
these data range from about 17 to 20 percent of the vacuum pressure coefficient. 

Shown in figure 7 are the variations of attainable thrust with normal Mach 
number for a series of airfoil sections. In spite of some inconsistencies, 
there is clearly a trend of reduced leading-edge thrust coefficients and reduced 
limiting pressure coefficients required for correlation as the Mach number 
increases. This study revealed a need for better and more complete two- 
dimensional airfoil data for symmetrical sections. In particular, a greater 
range of angle of attack and a greater range of Reynolds numbers would be 
desirable. 

A summary of the data used in defining the limiting pressure coefficient is 
shown in figure 8. The ratio of the limiting pressure to the vacuum pressure is 
shown as a function of the normal Mach number. There is obviously a strong 
dependence of the limiting pressure coefficient on the normal Mach number. 
Although the absolute value of the pressure coefficient decreases with increas- 
ing Mach number, the fraction of the vacuum pressure increases to values 
approaching 1.0 for Mach numbers near 1.0. There is also a weak but definite 
tendency of the limiting pressure coefficient to increase with increasing normal 
Reynolds number. The curve fit shown on figure 8 is intended to cover the Mach 
number and Reynolds number trends of the data. The curve fit is defined by the 
equation 


“'P , 1 im 


-2 


R n x 10 -6 


_ fi (4-3M n ) 

Rj, x 10 6 + 10 n 


0.05+0.35 (1-M n ) 2 


( 10 ) 


This form was chosen to allow the limiting pressure to approach the vacuum pres- 
sure as the Reynolds number approached infinity and to allow the limiting pres- 
sure to approach zero as the Reynolds number approached zero. With this 
definition of the limiting pressure, the system for estimation of attainable 
thrust is complete. 


NOTES ON METHOD IMPLEMENTATION 

The system for estimating attainable leading-edge thrust described in the 
preceding section of this report is intended for use as a subroutine in lifting- 
surface programs which provide estimates of theoretical leading-edge thrust. 

The following discussion covers program additions required to supply information 
needed in implementation of the system, reviews the steps in estimation of 
attainable thrust, describes the determination of flat-wing lift and drag coef- 
ficients with attainable thrust taken into account, and outlines the extension 
of the method to wings with twist and camber. 
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Computer Program Information Required 


In order to implement the present method as a subroutine in existing com- 
puter programs , several requirements must be met. The first requirement is for 
additional input data describing the streamwise airfoil sections. Thickness 
ratio t/c, leading-edge radius to chord ratio r/c, and position of maximum 
thickness g must be described as functions of span position. It is assumed 
that the wing-chord variation with spanwise position is already suitably 
described. Leading- and trailing-edge sweep angles, also, are expected to be 
provided or readily extracted from planform data. In addition, a Reynolds num- 
ber based on the wing mean aerodynamic chord must be input. 


Programming of Attainable Thrust Estimate 

In this review of the steps to be taken in estimating the local attainable 
spanwise thrust distribution, the original numbering of the equations developed 
in the section "Method Development" is retained. For each of a large number of 
selected span stations, a normal Mach number, normal section parameters, and a 
normal thrust coefficient must be calculated. These are, respectively. 


M n = m cos A le 


( 1 ) 


2nc 

C -|^% — * il. " * * ' ■ 1 k 1 l ^ = ^ ^ ~ — - — ■ — ■ — “ — .. 1 — ■ — — — 

sin A le j(l - n) tan A^ e + g tan A te j + cos A le 

Jn = T 1 

c n c 2r| cos A le 


( 2 ) 

(3) 


r n _ £ 1 

c n c 2ri cos 2 A]_ e 


"t , n 



1 

2 A 

cos A le 


(4) 

( 5 ) 


In addition, a normal Reynolds number and a limiting pressure coefficient corre- 
sponding to that Reynolds number and the normal Mach number are required: 


R n = R -f cos A le 


-2 


g? , 1 im 2 

Y M n 


R n x 10 -6 


R n x io “ 6 + io (4 3Mr) 


0. 05+0. 35 (1-M n ) 2 


( 6 ) 


( 10 ) 
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This then permits the calculation of an equivalent Mach number M 0 as 



(8) 


Finally, the thrust factor and the attainable thrust coefficient may be 
determined as 



— I 

M 0 ' 4 " 

2(l - M e 2) 

c n 


M e 

_ c t,n y/ 1 “ M n 2 _ 


but not greater than 1.0 


(9) 


The limitation of Krp to values no greater than 1.0 permits attainable thrust 
to equal but not exceed the theoretical thrust. 


Lift and Drag for a Flat Wing 

The integration of leading-edge thrust forces to obtain wing lift and drag 
coefficients may be handled so as to be compatible with the Polhamus leading- 
edge suction analogy (ref. 8). As shown in figure 9, the leading-edge suction 
vector is c^-/cos Aq e . In the Polhamus analogy which assumes that no leading- 
edge thrust is developed, this vector is assumed to rotate to a position per- 
pendicular to the wing surface and contribute to normal force rather than axial 
force. The rotated vector represents the effect of a fully developed separated 
leading-edge vortex system. In the present method where leading-edge thrust is 
partially developed, it would seem logical to consider a partial rotation of the 
vector. In this analysis, the vector is assumed to rotate only enough to give 
a chord-plane component equal to the predicted attainable suction, c£/cos A^ e . 
The component of the rotated vector perpendicular to the chord plane is assumed 
to give an incremental normal force associated with a partially developed sepa- 
rated leading-edge vortex system. In equation form. 


AC A = 


-c 


* 

t 


AC N " cos A 


sin cos 


-1 


le 


The incremental normal force vector is shown in figure 9 acting at the wing 
leading edge. Actually, the effect of the separated leading-edge vortex would 
be felt at some distance behind the leading edge. For a flat wing, the incre- 
mental normal force is not sensitive to this location, but the pitching moment 
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is. In this analysis, no attempt was made to define an effective location of 
the vector representing the separated vortex; it was assumed to act at the wing 
leading edge. 


Extension to Wings With Twist and Camber 


The present method for estimation of attainable leading-edge thrust has 
been developed for flat wings with symmetrical sections. However, the method 
should be adaptable to wings with limited twist and camber when it is coupled 
with lifting-surface programs capable of providing accurate theoretical leading- 
edge thrust distributions. Figure 10 will help to illustrate this application. 
Since the airfoil profile in the immediate vicinity of the leading edge has a 
dominant influence on the thrust characteristics, it should be possible to 
analyze the attainable thrust by performing calculations for a comparable sym- 
metric wing section. This section would have a plane of symmetry which is 
tangent to the mean camber surface of the nonsymmetrical section at some point 
(perhaps the center of the leading-edge radius) near the leading edge. The 
superimposed symmetrical section could be assumed to have the same thickness 
ratio, leading-edge radius, and location of maximum thickness as the cambered 
section. The thrust vector would be assumed to act at an angle 6 with respect 
to the wing-chord plane defined by the tangent to the camber surface. If the 
camber angle 6 is small, the incremental axial and normal force coefficients 
can be defined as 


where 


/-i t 1 

Ac A = “ c t cos 1 cos — - 6 ny 

1 . ( -1 c t fi \ 

t Sin ^ OS 1 ^7 - 6 nJ 


ACn — c 


tan 1 


/ tan 6 \ 

y os A le / 


These incremental section force coefficients could then be used in a manner 
similar to those for the flat wing to provide estimates of twisted- and cambered- 
wing aerodynamic performance for realistically attainable levels of leading-edge 
thrust . 

As for the flat wing, the incremental force normal to the lifting surface 
associated with a detached vortex system is assumed to act at the wing leading 
edge; whereas, it may actually be felt well aft of the leading edge. For a wing 
with twist and camber, not only will the wing normal force and pitching moment 
vary with the assumed vector location, but an additional incremental axial force 
will also be dependent on this location- Thus, the limited camber surface 
assumption (small 6) is quite restrictive. It should also be noted that care 
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must be taken to insure that the vector rotation for a twisted and cambered 
wing conforms to the direction of the local lifting forces at the wing leading 
edge. For twisted and cambered wings, there is clearly a need for a study of 
the proper handling in prediction techniques of the effects of partial thrust 
and the partial development of a separated vortex system. 


COMPARISONS WITH EXPERIMENTAL DATA 

The present method for estimation of attainable thrust was tested by 
application to wing-body configurations with flat wings and symmetrical sections 
for which experimental data (ref. 9) were available. Theoretical leading-edge 
thrust distributions were evaluated for subsonic and supersonic speeds by means 
of computer programs described in references 1 and 10, respectively. The attain- 
able thrust estimate was then obtained by input of these distributions into a 
special, separate computer program which implemented the computational process 
described in the section "Notes on Method Implementation." In computing the 
lift curve slope and the theoretical thrust distribution, the complete wing 
planform through the body was employed. This was done so that body carry-over 
lift would not be neglected in the calculation of total lift and in the deter- 
mination of leading-edge thrust. However, after theoretical leading-edge thrust 
distributions were obtained, that portion of the thrust inboard of the wing-body 
juncture was ignored. Theoretical lift and drag coefficients were formulated so 
as to be compatible with the Polhamus leading-edge suction analogy as described 
in the section "Notes on Method Implementation." 

The comparisons of theory and experiment shown in figures 11 to 16 are all 
of the same form. Axial force coefficient, angle of attack, drag coefficient, 
and lift-drag ratio are given as functions of the lift coefficient. The axial 
force is obviously the most sensitive indicator of the presence of leading-edge 
thrust and should be the primary gage of the success or failure of the estima- 
tion method. The lift-drag ratio is the single most important measure of the 
wing-body aerodynamic efficiency, in which leading-edge thrust plays a major 
role, and also is an important test for the system. In observing the lift 
coefficient-drag coefficient polar, the transition from full or near full 
thrust at low lift coefficients to smaller values of thrust at high lift coeffi- 
cients is of interest. The lift curve slope is defined by the basic lifting 
surface program, and the degree of leading-edge thrust has only a small influ- 
ence. For these correlations, no attempt was made to calculate theoretical 
zero-lift drag; instead, the experimental drag was used. To help assess the 
importance of leading-edge thrust, theoretical data are shown for the limiting 
conditions of zero and full theoretical thrust. 

In figure 11, results for the present method are compared with experimen- 
tal data for an unswept wing at a Mach number of 0.6. If the method has been 
properly formulated, it should work for this situation where the two-dimensional 
flow on which the system is based is a predominant factor. In spite of the 
rather thin wing section (3.0-percent thickness to chord ratio) and the very 
small leading-edge radius (0.045 percent of the chord), an appreciable fraction 
of the theoretical thrust is actually developed as the method predicts. 
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The remainder of the comparison figures are presented as pairs or as 
series in which one parameter varies as the others remain essentially constant. 

In part (a) of figure 12 , data are shown for a delta wing with a leading- 
edge sweep of 45.0° while in part (b) data are shown for a sweep angle of 63.4°. 
In both instances, the new method gives a good estimate of the measured aerody- 
namic characteristics at a Mach number of 0.25 for the 5-percent- thick wing. 

For the 63.4° swept wing, the prediction of lift and drag is seen to be reason- 
ably good in spite of some error in the axial force. 

The three parts of figure 13 show data for a 63.4° swept delta wing with 
thickness ratios of 3, 5, and 8 percent at a Mach number of 0.24 and a Reynolds 
number of about 5 x 10^. Except for the 8-percent-thick wing, there is no 
appreciable discrepancy between the experimental data and the estimate by the 
present method. For this sweep angle and this wing section, the normal airfoil 
thickness ratio is about 0.3; whereas, the largest thickness used in derivation 
of the thrust factors was 0.18. However, there is a reasonably good prediction 
of the lift-drag ratio in spite of the discrepancy between predicted and mea- 
sured thrust at the higher lift coefficients. 

The variation of attainable thrust with Mach number may be observed in the 
four parts of figure 14. The wing is a 63.4° swept delta with a 5-percent- thick 
section. At a Mach number of 0.24, the data are little different from that seen 
previously. For a Mach number of 0.6, the present method predicts a decrease 
in attainable thrust (compared to M = 0.24) which is not shown by the experi- 
mental data. Nevertheless, the lift and drag predictions are still good. At a 
Mach number of 0.9, the experimental data and the prediction show a decrease in 
attainable thrust, but this does not result in a decrease in the lift-drag ratio 
because of the higher lift curve slope. Because of the difficulty of predicting 
transonic flows, the good correlation shown here may be somewhat fortuitous. At 
the supersonic Mach number of 1.3, there is still evidence of some degree of 
leading-edge thrust. Although the thrust characteristics are well predicted, 
there is an appreciable difference in the lift-drag ratios. The discrepancy 
between the theoretical and measured lift at 0° angle of attack accounts for 
most of this difference. 

The variation of attainable thrust with Reynolds number may be explored 
with the aid of figure 15. The range of Reynolds numbers covered is not as 
large as desired; nevertheless, Reynolds number trends are evident. Thrust 
levels and trends indicated by the predicted axial force coefficient agree 
closely with the experimental data. There are, however, some discrepancies in 
the lift-drag ratio comparisons. There is some evidence of partially laminar 
flow and reduced skin friction for data at the lower Reynolds number, particu- 
larly at the lower lift coefficients. This may explain the behavior of the 
experimental data. 

In a further attempt to define limits of applicability, unpublished data 
obtained at three Mach numbers for a thick wing (NACA 0008 section) with a 
cranked leading edge are shown in figure 16. The data were obtained as part of 
the SHIPS (Subsonic/Hypersonic Irregular Planform Study) Program conducted by the 
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National Aeronautics and Space Administration. For the inboard portion of this 
wing, the normal thickness ratio is 0.5, a value much larger than any considered 
in derivation of the thrust factor. Since this configuration employed a non- 
symmetric fuselage, the angle of attack for zero lift as well as the zero lift 
drag used for the theory were taken from the experimental data. At the Mach 
number of 0.6, there is some overestimation of the thrust at the larger lift 
coefficients, but the lift-drag ratio estimate is still good. The prediction of 
thrust characteristics for a Mach number of 0.9 is rather poor, but for a Mach 
number of 1.2, the prediction by the present method agrees well with the experi- 
mental results. 


DESIGN APPLICATIONS 

Because the present method of predicting attainable leading-edge thrust 
takes into account the spanwise variation of airfoil section characteristics, 
an opportunity is afforded for design by iteration to maximize the attainable 
thrust and the attendant performance benefits. A simplified example of such use 
of the method is given in figure 17. A baseline wing-body configuration having 
a delta wing with aspect ratio of 4 was considered. This wing employed an 
NACA 0005-63 section throughout the span. The design conditions chosen were a 
Mach number of 0.6 and a lift coefficient of 0.26 at a wind-tunnel Reynolds 
number of 1.5 x 10 6 . As shown in the spanwise thrust distribution for the 
baseline configuration at the left of figure 17, that configuration is esti- 
mated to achieve about two- thirds of the theoretical thrust at design conditions. 
A wing such as this could not be expected to develop any appreciable thrust at 
the outboard span stations because of limitations imposed by the sharp tip where 
the chord and absolute thickness both go to zero. The dashed line shown in fig- 
ure 17 indicates the beginning of thrust limitation (this curve can be found 
from eq. (9) by setting K T = 1.0 and. solving for c t , n ) • Levels of theoreti- 
cal thrust anywhere below this line, where K<p is uniformly equal to 1.0, are 
estimated to be fully realizable. Above this line, K T is less than 1.0 and 
the theoretical thrust will not be fully developed. 

Redesign of the wing planform in an attempt to achieve a greater leading- 
edge thrust first involved removal of the limitations imposed by the zero tip 
chord, as indicated by the beginning of limitation line. When the planform was 
changed to that shown in the center of figure 17, the theoretical thrust level 
was reduced (only because the angle of attack corresponding to the design lift 
coefficient was reduced) , but the fraction of attainable thrust was increased 
to about 0.78 which resulted in an increase in attainable thrust. The increase 
in attainable thrust coupled with the improved lift curve slope is estimated to 
increase the lift-drag ratio for the design condition from about 15.0 to 16.7 
provided C d ^q remains unchanged. 

A further alteration in planform as shown at the right of figure 17 failed 
to produce an improvement over the trapezoidal planform. This does not mean, 
however, that a planform with curved leading edges producing a further improve- 
ment in aerodynamic efficiency cannot be found. 

Further improvement in the aerodynamic performance of the trapezoidal wing 
by means of a change in the spanwise thickness distribution is illustrated in 


16 



figure 18. Since theoretical thrust for the design condition is indicated to a 
span position of about 0.3Y/(b/2), thickness changes are required only out- 
board of this station. A linear increase in thickness ratio from 0.05 at the 
0.3y/(b/2) station to 0.08 at the tip station is estimated to increase the 
attainable thrust to about 96 percent of the full theoretical thrust. If the 
zero-lift drag is unaltered by this change, the lift-drag ratio for the design 
condition is estimated to be increased to about 19. 

Application of the present method to the design of wings for supersonic 
cruise vehicles may also be possible. Normally thrust considerations are 
ignored in supersonic aerodynamic theory, and wing lifting efficiency is opti- 
mized through use of twist and camber alone (e.g. , ref. 10) . The resultant 
wing camber surfaces, however, may be too severe for incorporation into prac- 
tical airplane designs. The large root chord angles and the resultant large 
cabin floor angles are particularly troublesome. If design by iteration using 
the present method could result in attainment of near theoretical leading-edge 
thrust over even a limited range of angle of attack or lift coefficient, the 
wing design lift coefficient and the resultant camber surface severity could be 
reduced accordingly with little or no loss in aerodynamic efficiency. 


CONCLUDING REMARKS 

A study has been made of the factors which place limits on the theoretical 
leading-edge thrust, and an empirical method for estimating attainable thrust 
has been developed. The method is based on the use of simple sweep theory to 
permit a two-dimensional analysis, the use of theoretical airfoil programs to 
define thrust dependence on local geometric characteristics, and the examination 
of experimental two-dimensional airfoil data to define limitations imposed by 
local Mach numbers and Reynolds numbers. The applicability of the method was 
demonstrated by comparisons of theoretical and experimental aerodynamic charac- 
teristics for a series of wing-body configurations. Generally, good predictions 
of the attainable thrust and its influence on lift and drag characteristics were 
obtained over a range of Mach numbers from 0.24 to 1.3. 

The method is compatible with the Polhamus leading-edge suction analogy for 
fully detached vortex flow, when the analogy is considered to represent the 
limiting condition of a gradual rotation of the suction vector that occurs as 
leading-edge thrust is lost. An additional advantage of the method is the 
possibility it presents for designing wings in an iterative manner to maximize 
thrust and the attendant performance benefits. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 24, 1979 
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APPENDIX 


DERIVATION OF NORMAL AIRFOIL SECTION PARAMETERS 

Derivation of normal airfoil sections and the appropriate flow parameters 
is based on simple sweep theory and on the use of a superimposed arrow-wing 
planform as shown in figure 1. This phantom wing has the same sweep angles for 
the leading edge, the trailing edge, and the maximum thickness line as does the 
actual wing at a given span station. 

The normal Mach number is simply the component of the free- stream Mach 
number in a direction perpendicular to the wing leading edge, as shown in 
sketch (Al) . 



M n = M COS A^ e (1) 

The chord of the normal section is defined to place the maximum thickness 
at midchord. The sweep angle of the maximum thickness line may be obtained by 
consideration of the superimposed arrow-wing planform shown in sketch (A2) . 
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The relationship between the normal section chord and the streamwise section 
chord may be established from the triangle shown in sketch (A3) . 
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APPENDIX 


From sketch (A4) , 



an expression for the maximum thickness of the normal section can be found: 
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APPENDIX 


A relationship between the leading-edge radius of the normal section and 
the leading-edge radius of the streamwise section can be derived by considering 
the streamwise airfoil section to be represented by the equation 

z - k iv pr + ^ k i+1 *< i+1 >/ 2 

in which only the first term contributes to the leading-edge radius, r = k -^f 2. 
In the normal plane. 


L n , | rjc , . , 

z n = — k 1 yx n — - + hxgher order terms 


k-,. — \/xT + higher order terms 

T l\0.5c n V n 


with a leading-edge radius of 


Then, 



(4) 
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APPENDIX 


The section thrust coefficient may be defined as the thrusting force per 
unit dynamic pressure, per unit chord, and per unit spanwise distance. Thus 
for the streamwise section, 

c = — -L 
t dy qc 


Sketch (A5) will aid in the definition of the normal section thrust coefficient. 



Sketch (A5 ) 


The thrust vector in the normal plane is 

dt 


dt n = 


cos 


A 


le 


The dynamic pressure in the normal plane is 


- - 


lr. - 'Am; * <5 cos2 S le 
and the incremental distance in the spanwise direction for the normal section is 
dy 


ds = 


v le 
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The normal section thrust coefficient then becomes 
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dt n 1 
ds Sn c n 


dt cos A le 1 

cos A le dy q cos 2 AieCn 

dt c_ 1 

d y <I C c n cos 2 A le 

c 1 

c t c 2 a 

c n cos z A le 


( 5 ) 


The normal section Reynolds 
number because of changes in the 
Thus , 


number differs from the streamwise 
velocity or Mach number and in the 


Reynolds 

chord. 


Rn = R T 


c n M n 


a m 


R IT- cos A 
c 


le 


( 6 ) 


23 



REFERENCES 


1. Lamar, John E.? and Gloss, Blair B. : Subsonic Aerodynamic Characteristics 

of Interacting Lifting Surfaces With Separated Flow Around Sharp Edges 
Predicted by a Vortex-Lattice Method. NASA TN D-7921, 1975. 

2. Jones, Robert T. : Estimated Lift-Drag Ratios at Supersonic Speed. NACA 

TN 1350, 1947. 

3. Carlson, Harry W. ; and Mack, Robert J. : Estimation of Leading-Edge Thrust 

for Supersonic Wings of Arbitrary Planform. NASA TP-1270, 1978. 

4. Henderson, William P. : Studies of Various Factors Affecting Drag Due to 

Lift at Subsonic Speeds. NASA TN D-3584, 1966. 

5. Bauer, Frances; Garabedian, Paul; Korn, David; and Jameson, Antony: Super- 

critical Wing Sections II. Volume 108 of Lecture Notes in Economics and 
Mathematical Systems, Springer-Verlag , c.1975. 

6. Daley, Bernard N. ; and Dick, Richard S-: Effect of Thickness, Camber, and 

Thickness Distribution on Airfoil Characteristics at Mach Numbers up to 
1.0. NACA TN 3607, 1956. (Supersedes NACA RM L52G31a.) 

7. Loftin, Laurence K. , Jr.: Aerodynamic Characteristics of the NACA 64-010 

and 0010-1.10 40/1.051 Airfoil Sections at Mach Numbers From 0.30 to 0.85 
and Reynolds Numbers From 4.0 x 10 6 to 8.0 x 10 6 . NACA TN 3244, 1954. 

8. Polhamus, Edward C. : Predictions of Vortex-Lift Characteristics by a 

Leading-Edge Suction Analogy. J. Aircr. , vol. 8, no. 4, Apr. 1971, 
pp. 193-199. 

9. Hall, Charles F. : Lift, Drag, and Pitching Moment of Low-Aspect-Ratio Wings 

at Subsonic and Supersonic Speeds. NACA RM A53A30, 1953. 

10. Carlson, Harry W. ; and Miller, David S.: Numerical Methods for the Design 

and Analysis of Wings at Supersonic Speeds. NASA TN D-7713, 1974. 


24 



Superimposed arrow 



Figure 1.- Relationship between streamwise and normal wing sections 
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Figure 2.- Airfoils employed in analytic study. 
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Figure 5.- Characteristics of two-dimensional airfoil for different 
thickness ratios; M = 0.3; R = 0.67 x 10 6 . 



O Experiment (ref. 5), large duct 
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Figure 6.- Characteristics of two-dimensional airfoil for different 
locations of maximum thickness; M = 0.3; R = 0.67 x 10^. 
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(a) NACA 64A004 airfoil. 

Figure 7.- Characteristics of two-dimensional airfoil for different 
Mach numbers and Reynolds numbers . 
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Figure 7.- Continued. 
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Figure 7.- Continued. 


by iteration 

_ M =0.6 

R = 1.20 X 10 c 



^ M = 0.9 

R = 1.54 x 10* 


'p,lim 

J L 


= -1.5 

\ 


2 4 6 

a, deg 


35 





I 


O Experiment (ref. 6) 
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Figure 16.- Comparison of theory and experiment for three Mach numbers; 
thick wing with cranked leading edge. 
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Figure 17. Example of use of present method for wing planform design. 




Figure 18.- Example of use of present method for airfoil section selection 
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